Saturation in heteronuclear photoassociation of 6 Li 7 Li 



(N 
O 

o 

> 
O 

(N 



U. Schloder, C. Silber, T. Deuschle, and C. Zimmermann 

Phystkalisches Institut, Eberhard Karls Universitat, 
Auf der Morgenstelle 14, 72076 Tubingen, Germany 
(Dated: February 9, 2008) 

We report heteronuclear photoassociation spectroscopy in a mixture of magneto-optically trapped 
6 Li and 7 Li. The laser-induced decrease in the 7 Li steady-state particle number, only appearing 
in the presence of 6 Li, gives clear evidence of photoassociation to form 6 Li 7 Li. Hyperfine resolved 
spectra of the vibrational level w = 83 of the singlet state A x Ej have been taken up to intensities of 
1000 W/cm 2 . The absolute resonance frequencies and the rotational constant have been measured. 
Saturation of the photoassociation rate has been observed for two hyperfine transitions, which can be 
shown to be due to saturation of the rate coefficient near the unitarity limit. Saturation intensities 
on the order of 40 W/cm 2 can be determined. 

PACS numbers: 32.80.Pj, 33.80.Ps, 33.70.Jg, 03.75.Fi 
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Photoassociation of cold atoms has been established as 
a versatile tool for the formation of cold molecules. Start- 
ing with atoms either from a magneto-optical trap or a 
Bose-Einstein condensate this technique has been shown 
to allow for efficient production of various cold alkali 
dimers (D, 0, |l 0]. I n single-photon photoassociation a 
pair of free atoms is optically excited into a rovibrational 
level of the electronically excited molecule. By sponta- 
neous decay, either a pair of free atoms is formed, or, 
with small fractional probability, a translationally cold 
ground state molecule. In extension to this intrinsically 
incoherent process, also coherent processes, such as two- 
photon Raman- type transitions {§ |[ @], have been 
investigated. In the degenerate regime, such methods 
should allow for coherent coupling between free atom- 
pair states and specific molecular ground states, similar 
to the recently observed atom-molecule coupling due to 
a Feshbach resonance ||. It is expected that quantum 
statistics will affect the initial atomic states as well as 
the molecular product, leading to what may be phrased 
as 'superchemistry' [Q. In this context, mixtures between 
bosonic and fermionic gases are particularly interesting, 
as they offer not only the possibility of molecular Bose 
gases, but also of molecular Fermi gases. 

A promising candidate for such experiments is lithium. 
It provides a stable fermionic and bosonic isotope, it 
can be routinely cooled to quantum degeneracy, and 
mixtures between Bose- and Fermi-gases have already 
been observed and studied |ll], pd[ . There is also ex- 
tensive experience with photoassociation of homonuclear 
lithium dimers in the nondegenerate and in the de- 
generate regime |l2|, However, the formation of 
fermionic molecules requires two atomic species with dif- 
ferent quantum statistics. In general, such heteronu- 
clear photoassociation bears intrinsic difficulties |I3 and, 
although being very desirable as a universal tool Jl5| , 
it has not yet been implemented successfully in binary 



mixtures. For studying heteronuclear photoassociation, 
lithium is again exceptionally suitable. The isotopic shift 
of only 10 GHz allows for resonant dipole-dipole interac- 
tion in the excited molecular state already at large in- 
ternuclear separations, which is a crucial prerequisite for 
efficient photoassociation fig] . 

Even though the theory of photoassociation at high 
laser intensities is not yet fully understood Jl7], |lj] , it is 
clear that for coherent coupling, optical excitation rates 
are required that are fast as compared to the other rel- 
evant time scales of the system (i.e. collision time and 
decoherence time due to spontaneous decay of the in- 
volved excited states). One indication for entering the 
regime of sufficiently strong excitation rates is the occur- 
rence of saturation in the single-photon photoassociation 
signal. Saturation studies have been reported for several 
one-species experiments ^0). However, a clear satu- 
ration of the photoassociation rate has not been observed 
yet. 

In this paper we report on heteronuclear photoasso- 
ciation in a mixture of magneto-optically trapped 6 Li 
and 7 Li. We investigate photoassociation of 6 Li 7 Li by 
monitoring the 6 Li-induced decrease in the 7 Li steady- 
state particle number as a function of the photoassoci- 
ation laser frequency and intensity. The number of 6 Li 
atoms exceeds that of 7 Li by a factor of 10 such that 
the photoassociation losses are significant only for 7 Li 
while the number of 6 Li atoms remains almost unaffected. 
This setup is possible only for heteronuclear photoassoci- 
ation and substantially facilitates the quantitative anal- 
ysis. Strong saturation of the photoassociation rate is 
observed for two of four studied hyperfine lines of the 
chosen singlet transition t? = 83. We derive values for the 
saturation intensities on the order of 40 W/cm 2 which are 
sufficiently low to encourage future experiments aiming 
at coherent control in the degenerate regime. 

Our experimental setup has been described in detail 
in Q. The combined magneto-optical trap (MOT) is 
operated with two independent diode-laser systems and 
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loaded from a single Zeeman-slowed atomic beam. The 
particle numbers are monitored by two absorption lasers, 
tuned into resonance of the respective cooling transition, 
by means of a lock-in technique. The photoassociation 
light is provided by a tunable dye laser, which is fo- 
cused to a 1/e 2 diameter of 300 /im, resulting in peak 
intensities up to 1000 W/cm 2 . Its frequency is measured 
with a combination of a commercial wavemeter, a Fabry- 
Perot etalon and a doppler-free iodine fluorescence spec- 
troscopy ]2lj . If the laser light is resonant with a molecu- 
lar transition, losses in the two-species MOT are induced 
due to spontaneous decay of the excited molecule into a 
pair of free atoms or a ground state molecule. To simplify 
the quantitative analysis, we operate the MOTs such that 
the particle number of the 6 Li MOT is large compared to 
that of the 7 Li MOT. Therefore, losses in the 6 Li particle 
number due to heteronuclear photoassociation are small 
and the 6 Li MOT acts as a temporally constant back- 
ground gas for the 7 Li MOT. Moreover, quadratic losses 
in the 7 Li particle number are negligible. In that special 
case the complicated solution of the coupled rate equa- 
tions for a two-species system is reduced to the following 
simple form for the 7 Li steady-state particle number N?: 



Here, L-j is the 7 Li MOT loading rate, and oc? is the loss 
rate mainly due to collisions with hot atoms from the 
beam source. The photoassociation rate /3pa = Kn^wr, 
with the photoassociation rate coefficient K, the 6 Li peak 
density uq, the factor w considering the relative atomic 
hyperfine ground-state populations, and the factor r 
accounting for the spatial overlap between the atomic 
clouds and the laser beam. Spectra are taken by moni- 
toring the 7 Li steady-state atom number as a function of 
the photoassociation laser frequency. With 1.25MHz/s, 
the scan rate was chosen slow as compared to the load- 
ing time of the MOTs. The positions of the heteronu- 
clear resonances can be predicted theoretically within 
1 GHz, by applying the method of mass-reduced quan- 
tum numbers |22[ to the well-known data for homonu- 
clear lithium [{l2f . To identify heteronuclear resonances 
we have recorded the 7 Li particle number with and with- 
out the other isotope present. The additional resonances 
give clear evidence of 6 Li 7 Li photoassociation. 

We have concentrated on the transition into the vi- 
brational level v = 83 of the singlet series A : E+ near a 
detuning of 177 GHz to the red of the Dl line of 6 Li. Its 
outer turning point is placed at an internuclear separa- 
tion of about 74 a . At this detuning, the singlet series 
of 6 Li 7 Li is sufficiently separated from the heteronuclear 
triplet resonances as well as from the homonuclear 6 Li 6 Li 
and Li Li transitions, while at the same time the pho- 
toassociation rates are still large. A high-resolution pho- 
toassociation spectrum for an intensity of 420 W/cm 2 is 
shown in Fig. El The four strongest lines ((a)-(d)) can be 
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FIG. 1: Heteronuclear photoassociation spectrum of the sin- 
glet v = 83 level of 6 Li 7 Li. The frequency scale is rela- 
tive to the hyperfine center of gravity, which amounts to 
14897.4015(7) cm -1 . The inset shows the assignments of the 
different hyperfine transitions. 



assigned to transitions from the quartet of hyperfine en- 
trance channels into the rotational level N = 1 of the ex- 
cited state. This hyperfine structure reflects the different 
combinations of atomic hyperfine ground states (hyper- 
fine splitting: 228 MHz for 6 Li and 804 MHz for 7 Li) and 
is a further clear indication of the heteronuclear struc- 
ture of the molecule. The absolute transition frequency, 
adjusted to the hyperfine center of gravity, amounts to 
14897.4015(7) cm -1 . Several lines to the excited-state 
levels N = (resonances (a')-(d')) and N = 2 (resonance 
(a")) are also visible. From more sensitive spectra taken 
with a less focused laser beam (1.0 mm diameter), the ro- 
tational constant B/h can be determined as 301(2) MHz. 
The transitions into the N = 1 level can be assigned to s- 
wave collisions and the transitions into the levels N = 0, 2 
to p-wave collisions. The contributing partial waves are 
determined in their parity by the selection rule on the 
change of angular momentum and in their number by 
the corresponding centrifugal barriers. For d-wave col- 
lisions the barrier amounts to 38 mK at an internuclear 
separation of 54 ao, which strongly suppresses contribu- 
tions from d-waves and above, even as shape resonances. 

In the following we discuss the intensity dependence 
of the photoassociation rate. We took several series of 
spectra up to intensities of 1000 W/cm 2 and find, that 
the signal saturates. For all hyperfine transitions, we ob- 
serve a remaining fraction of 7 Li atoms, that cannot be 
depleted, even in the limit of high laser intensities. Thus, 
particle losses due to photoassociation must be limited at 
high intensities. This limitation is to be explained by the 
intensity dependence of the photoassociation rate. Inter- 
nal MOT processes, such as optical pumping, thermal- 
ization or diffusion, can be excluded as the origin of the 
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FIG. 2: Spectra for the f 6=3/2 + f 7=2 transition (resonance 
(a)) at two different intensities. The solid line (lower panel) 
is a fit according to Eq. ^, the corresponding residuals are 
plotted in the upper panel. The frequency scale is relative to 
the central frequency. 



observed behaviour by a simple time scale argument. In 
order to lead to a limited depletion, the losses due to pho- 
toassociation must occur on the same time scale as the 
loading time of the MOT, which is about 10 s. Internal 
MOT processes happen on a much faster time scale and 
therefore cannot limit the losses. 

For a quantitative analysis, one has to take into ac- 
count that the experimentally observed steady-state par- 
ticle number Nj depends on the photoassociation rate in 
a nonlinear way (see Eq.|). We assume a Lorentzian fre- 
quency dependence of the photoassociation rate /3pa = 
a/((/ - fres) 2 + 7 2 / 4 ), with the central frequency f res , 
the linewidth 7 and a proportionality factor a. Inserting 
this expression in equation |l| yields: 



iV 7 



JV7,0(1 - Cf) 



l + a/(a r ((/-/ res ) 2 +7 2 /4)) 



(2) 



We account for slow drifts in the loading rate during 
the frequency scan and for partial overlap of the dif- 
ferent lines by assuming a linear frequency dependence 
oc (1 — cf) of the steady-state particle number without 
photoassociation laser present 2VV 0- For the atomic-beam 
limited loss rate aj we use the measured value of 0.11 s _1 . 
Up to intensities of 450 W/cm 2 for transitions (a) and (b) 
and of 600 W/cm 2 for transitions (c) and (d) the different 
hyperfine lines do not show any significant asymmetries 
and partial overlap of the different lines remains small, 
so that the resulting quasi-Lorentzian function represents 
the data sufficiently well. This is exemplarily shown in 
Fig. U for transition (a) at two intensities. Obviously, 
asymmetries due to thermal averaging over the differ- 
ent collision energies |2J| are negligible, since the natural 
linewidth of the excited level (12 MHz, ||]) plus the 
unresolved hyperfine-structure of the excited state [Q is 
larger than the thermal broadening of about 10 MHz for 
the MOT temperature of 0.5 mK p7|. 
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FIG. 3: Photoassociation rate as a function of the laser peak 
intensity for the four different hyperfine transitions. The open 
circles represent the uncorrected data. 



For the determination of the photoassociation rate one 
has to realize that the shape of the resonances can be 
influenced by various systematic effects. Due to the fo- 
cussing of the photoassociation laser and due to inho- 
mogenities in the beam profile, the light shift jl£| varies 
spatially, which leads to inhomogeneous broadening. The 
different light shifts for the different hyperfine levels of 
the excited state contribute similarily. These broadening 
mechanisms lead to asymmetries at very high intensities 
and flatten the approximately Lorentzian lines at moder- 
ate laser intensities. However, the area of the frequency- 
integrated photoassociation rate is preserved, which has 
been checked numerically. Therefore, we correct the max- 
imum photoassociation rate 4a/7 2 , which we get from 
the fit, by multiplying with the ratio of the broadened 
linewidth 7 and the linewidth in the limit for zero in- 
tensity, 77=0, which gives /?pa,o = 4a/77/ =0 p8j. This 
corrected photoassociation rate /?pa,o is proportional to 
the area of the rate and corresponds to the maximum rate 
in absence of any broadening mechanisms, both broaden- 
ing due to systematic effects and power broadening due 
to saturation. Extrapolation of the broadened linewidths 
towards zero intensity leads to values for 77=0 of 51, 45, 
33, 40 MHz for resonances (a)-(d) respectively. They in- 
clude the natural linewidth, the unresolved structure of 
the different involved excited hyperfine states and ther- 
mal broadening. In Fig. || the photoassociation rate for 
the four different hyperfine transitions is presented as a 
function of the laser intensity. The data points are mean 
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values of up to 4 spectra. The error bars indicate the 
statistical error. While the rates for the transitions (c) 
and (d) do not saturate in the given intensity range (the 
uncorrected rates 4a/7 2 , which seem to saturate are also 
indicated), strong saturation is observed for the transi- 
tions (a) and (b). As the 6 Li density is nearly unaffected 
by the photoassociation, this saturation must be due to 
saturation of the rate coefficient K. 

For the interpretation we use a simplified model ac- 
cording to the close-coupled theory of Bohn and Juli- 
enne ft yields an analytic expression for the rate 
coefficient for a two-level system and a specific scatter- 
ing energy. To compare with the experimental corrected 
photoassociation rate, we have to adapt the theoretical 
expression by multiplying with the ratio of the power- 
broadened linewidth and the natural linewidth. For s- 
wave collisions this leads to the following intensity depen- 
dence for the corrected photoassociation rate coefficient: 



nv 



k 2 1+Isat/I 



(3) 



where v is the velocity and k is the wave number of the 
relative motion. I sat is the energy-dependent saturation 
intensity, for which the natural linewidth equals the in- 
duced linewidth. This intensity dependence also applies 
for the thermally-averaged corrected rate coefficient with 
approximately the same saturation intensity, as checked 
by simulations. A full theoretical analysis would include 
a population- weighted summation over the different tran- 
sitions from the magnetic ground-state sublevels to the 
various excited hyperfine states and has not been carried 
out yet. We thus assume the validity of the expression 
(Eq. ||) also for the full problem, and use it for a fit to the 
data of the saturating transitions. For the saturation in- 
tensities we obtain values of 28 W/cm 2 for transition (a) 
and of 54 W/cm 2 for transition (b). This is about 10 4 
times the atomic saturation intensity and in good agree- 
ment with an estimation for a two-level system, which 
gives 33W/cm 2 for an energy of E/k-g — 0.5mK [ p9| . 
For the corrected photoassociation rates, we receive max- 
imum values of 0.06s -1 for resonance (a) and O.I4s _1 
for resonance (b). Taking into account an estimated 6 Li 
peak intensity of 10 10 cm -3 , the overlap factor r = 0.22, 
the population factor w (0.42 and 0.2f for the transitions 
(a) and (b) respectively, statistical distribution assumed) , 
and the correction factor of 4 (see Eq. ||) , the saturated 
photoassociation rate coefficient K can be estimated to 
be (2-10~ n - 8-10" 11 ) cm 3 s _1 . This value is close to the 
unitarity limit, where the scattering probability \S\ 2 — > 
1. In this limit the rate coefficient is simply given by the 
product of the scattering cross-section ir/k 2 and the ve- 
locity v, resulting in 5-1 -10 cm 3 s _1 for the temperature 
considered. For a more sophisticated analysis, knowl- 
edge of the different ground-state sublevel populations, 
of the various Clebsch-Gordan coefficients and of laser- 
polarization effects would be required. 



In summary, we have observed heteronuclear photoas- 
sociation in 6 Li 7 Li. Studying the intensity dependence 
of the hyperfine resolved singlet spectrum, we have ob- 
served saturation of the rate coefficients for two of the 
four hyperfine lines and determined the corresponding 
saturation intensities. The observation of the unitarity 
limit has been possible, as it is orders of magnitude lower 
at MOT temperatures than at temperatures in the de- 
generate regime ]2Cj . Reaching the unitarity limit is an- 
other proof that already at MOT temperatures quantum 
mechanical properties become apparent. 
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